Interfacial Effects of Al- Termination on Spin Transport in Magnetic Tunnel Junctions 
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Experiments have shown that the tunneling current in a C0/AI2O3 magnetic tunneling junction 
(MTJ) is positively spin polarized, opposite to what is intuitively expected from standard tunneling 
theory which gives the spin polarization as exclusively dependent on the density of states (DOS) at 
Ef of the Co layers. Here we report theoretical results that give a positive tunneling spin polarization 
and tunneling magnetoresistance (TMR) that is in good agreement with experiments. From density 
functional theory (DFT) calculations, an Al-rich interface MTJ with atomic-level disorder is shown 
to have a positively polarized DOS near the interface. We also provide an atomic model calculation 
which gives insights into the source of the positive polarization. A layer and spin dependent effective 
mass model, using values extracted from the DFT results, is then used to calculate the tunneling 
current, which shows positive spin polarization. Finally, we calculate the TMR from the tunneling 
spin polarization which shows good agreement with experiments. 

PACS numbers: 



I. INTRODUCTION 

The discovery of large tunneling magnetoresistance 
(TMR) in magnetic tunneling junctions 1 (MTJs) has 
prompted great interest in developing these devices for 
applications such as magnetic random access memo- 
ries (MRAMs) and magnetic sensors (for reviews, see 
e.g. Refs. [2.3]). The TMR should, in an ideal system, 
be closely related to the spin polarization of the two fer- 
romagnetic layers sandwiching the insulator, and is typ- 
ically given by the Julliere model^ as 



TMR = 



Rap — Rp 
Ra~p 



2PiP 2 
1 + P1P2 



(1) 



Here Rap and Rp denote the resistance for the anti- 
parallel and parallel alignments of the ferromagnetic lay- 
ers, respectively, and Pi (i = 1 and 2) are the bulk mag- 
netic moments of the two layers. This model is based on 
the assumption that the spin polarization of the tunnel- 
ing current, which determines the TMR, is an intrinsic 
property of the spin polarization in the ferromagnetic lay- 
ers. Experimental results, however, do not often match 
this phenomenology. Most importantly, for the ferromag- 
netic 3d metals such as Co, which have a negative spin 
polarization at the Fermi level, the measured spin po- 
larization across an alumina, AI2O3, insulating barrier is 
positive. 5 . Experiments have also demonstrated a strong 
dependence of spin transport on the detailed structural 
and electronic nature of the insulating layer and its in- 
terface with the ferromagnetic layers^i 

Density functional theory (DFT) calculations of ex- 
terminated interfaces have shown a positive spin polar- 
ization for barrier distances beyond 10 A4>£. Additional 
DFT work have shown that additional oxygen atoms 
absorbed at the interface form strong Co-0 bonds and 
give rise to a positive spin interface band. This band 
was found to be dominant in the tunneling process, and 



gives rise to a positive spin polarization of the tunneling 
current. 10 Formation of an O-terminated interface with 
Co-0 bonds is obtainable by moderate to long oxidation 
times of the AI2O3 layer. An increased magnetization of 
the interface Co atoms was found to coincide with the 
optimum TMR, using X-ray magnetic circular dichroism 
measurements M- However, one recent experiment failed 
to detect any induced magnetic moment on the interface 
O-atomsjiSi which should also be present. Hence, details 
of the physical nature of O-terminated MTJs remains 
open. 

In this work we will instead focus on clean Al- 
tcrminatcd and disordered Al-rich interfaces which 
are experimentally accessible by not overoxidizing the 
junction^ Earlier DFT results on Al-terminated junc- 
tions have shown that the magnetic moment on the Co 
interface atoms is somewhat reduced by charge trans- 
fer and screening effects^ but Al-rich junctions have not 
yet been studied nor have the positive spin polarization 
been explained in any of these junctions. We show here 
that increasing Al content at the interface can cause dra- 
matic changes in the spin polarization of both the in- 
terface region and the tunneling current, even to the 
degree of changing the sign of the spin polarized tun- 
neling current as measured experimentally This explic- 
itly demonstrates that oxygen absorption is not the only 
means to achieve positive spin polarization as was re- 
cently suggested^ 

More specifically, we use DFT to model C0/AI2O3/C0 
magnetic tunnel junctions with varying Al-rich interfaces 
to determine atomic structure, band structure, and lo- 
cal density of states (LDOS). The DFT results explic- 
itly show the importance of the proximity effect on the 
LDOS of the interfacial layers. The LDOS near Ep of 
the interfacial Co layer leaks into the first Al layer of 
the insulating barrier, and gives rise to a small but finite 
DOS near Ep, hence metallizing the Al layer, as has been 
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seen experimentally^ The Al layer also picks up a small 
ferromagnetic moment due to this proximity effect. Sim- 
ilarly, the LDOS and the spin polarization of the nearest 
Co layer is reduced. Thus, the spin DOS changes as we 
cross the interface, and the spin polarization actually be- 
comes positive for Al-rich interfaces. 

We use a position-dependent effective mass model to 
explicitly calculate the spin polarized tunneling current 
of our junctions. The change in the effective mass of the 
different bands between the electrode bulk-like layers and 
the interfacial layers result in better matching between 
the spin- up bands, compared to the spin-down bands. 
There is then less reflection at the interface, and hence 
a larger transmission coefficient for the spin-up bands 
ultimately resulting in a positive spin polarized tunneling 
current in the Al-rich junctions. 

The paper is organized as follows. In Section [II] we re- 
view the details specific to the DFT calculations. We also 
present the different interface structures we have stud- 
ied. In Section IIIII we report the results for the clean 
Al-terminated interface as a baseline for the Al-rich in- 
terfaces. In Section IIVI we discuss interfacial scattering 
by excess Al atoms before we present DFT results for 
disordered, Al-rich interfaces in Section fVj Finally in 
Section IVTl we explicitly calculate the spin polarized tun- 
neling current and the TMR based on the DFT LDOS 
and band structure data. 



II. DFT COMPUTATIONAL DETAILS 

To model the C0/AI2O3/C0 magnetic tunnel junction 
we have created C0/AI2O3/C0 supercells with various 
Al-terminated interfaces. Thin film Co is predominantly 
found in the fee phase and experimentally it has been 
found that alumina grows on top of the (111) plane of 
fee Co.— We have chosen a 2x2 surface unit cell of (111) 
fee Co with the theoretical lattice parameter a = 3.38 A 
(a-exp = 3.55 A) as the base structure. In order to en- 
force bulk conditions in the inner cobalt layers we fix 
the atomic positions for the three innermost Co lay- 
ers. For the tunneling barrier we use one unit cell of 
a-Al 2 03 (corundum) with the [0001] orientation and Al- 
tcrmination to model the experimentally amorphous alu- 
mina. This structure contains four Al layers with three 
O layers in between them. Since alumina is grown on top 
of cobalt we adjust the lateral dimensions of alumina to 
that of the cobalt. This is the same crystallographic ori- 
entation as used in previous DFT studies^^ and gives 
an experimental lattice mismatch of only 6% between 
the cobalt and alumina. While most previous theoret- 
ical work have focused on oxygen rich interfaces there 
also exist experimental data showing Al-termination at 
the C0/AI2O3 interfaced We investigate both the clean, 
abrupt interface with Al-termination, see Fig. [TJ and two 
different cases of Al-rich interfaces. The first case has one 
mixed interface layer with a 1:1 Al to Co ratio. The sec- 
ond case has two mixed layers with a 3:1 and 1:1 Co to Al 



ratio, respectively, see Fig. [5l We will present the results 
for the clean interface system in Sec. IIIII and the results 
for the disordered interface system with two mixed layers 
in Sec. [V] We found the results for the system with one 
mixed layer to be an interpolation of the clean and two 
mixed layer interfaces. 

The band structure and LDOS were calculated using 
the first-principle density functional theory (DFT) pseu- 
dopotential method implemented in the Vienna Ab-initio 
Simulation Package (VASP),— We have used ultrasoft 
pseudopotentials and employed both the local density 
approximation (LDA) and the general gradient approxi- 
mation (GGA). In general, we expect the LDA calcula- 
tions to underestimate the lattice constant and the band 
gap but with more consistent errors than a GGA calcu- 
lation. The cut-off energy for the plane wave expansion 
was set to 29 Ry. After k-point convergence tests we 
chose a 9x9x1 T-centered k-point sampling. The atomic 
structure, except the fixed three innermost Co layers, 
was relaxed until the change in total energy between two 
ionic steps was less than 10~ 3 eV. The LDOS per layer 
was calculated as a spd site-projected LDOS using the 
monoatomic Wigner-Seitz radii and then summed over 
each layer, see Figs. [1] and [5] for typical layout of the 
layers. 




FIG. 1: Atomic structure of the clean Al-terminated interface 
junction. The unit vectors for the supercell as well as the layer 
identification tags are shown. Layers Co(l), Co(2), and Co(7) 
were kept fixed in order to impose bulk Co conditions. 
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III. CLEAN INTERFACE 

Figure Q] shows the supercell used for calculating the 
properties of a clean, Al-terminated interface junction. 
We performed both LDA and GGA calculation but found 
the LDA results to be more self-consistent. For reference 
we also used Co in the hep structure but could not detect 
any significant change from the fee structure, which was 
then used throughout the rest of the work. 

Figures [2][4] show the LDOS for the innermost Co layer, 
Co(l), the interface Co, Co(4), layer, as well as the first 
Al layer, Al(l). As expected, the LDOS for Co(l), which 
is furthest from the interface, is almost identical to that 
of bulk Co with an exchange splitting of « 1.5 eV. The 
major peaks in the spin-up d-electron LDOS is below 
Ep , indicating that the spin-up d-orbitals are filled as 
expected. Similarly, some of the spin-down d-electron 
peaks are above Ep, showing that some of the spin-down 
d-orbitals are unfilled; this gives rise to the net mag- 
netic moment. The s- and p-electrons also show some 
hybridization with the d-electrons, as evidenced by the 
slight increase in the spin-up s- and p-electron density 
at w 3 eV, and the spin-down s-electron and p-electron 
density at » 4.3 eV. 
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FIG. 2: LDOS for the innermost Co layer, Co(l), for the clean 
Al-terminated interface. The LDOS shows typical bulk-like 
character as expected. 



The spin polarization, <5 S = ^|| £ ,^j + ^|| £ ,^ , for the Co(l) 
layer is -0.64, but -0.43 for the Co(4) layer. For this clean 
interface system, this is due to the hybridization between 
the Co d-electrons and the Al s-p electrons at the inter- 
face. The trend of increasing positive spin polarization as 
the interface is approached is seen in all three systems, 
i.e. in both the clean interface and the two disordered 
interface cases which we have modeled. The detailed be- 
havior of the spin polarization is shown in Fig. 1 101 and 
is discussed further in Section. [V] 
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FIG. 3: LDOS for the interfacial Co layer, Co(4), for the clean 
Al-terminated interface. 

The Al(l) layer is the start of the insulating AI2O3 
barrier. However, due to hybridization of the Al s-p 
electrons with the d-electrons in the neighbouring Co(4) 
layer, there is a finite DOS at Ep. It is only the 0(1) 
layer next to Al(l) that shows a gap at Ep and therefore 
the actual width of the insulating barrier is decreased by 
one to two atomic layers. Calculation of the current- 
voltage curves using both cases show that a decrease 
of one atomic layer gives a better fit with experimen- 
tal results^ This means that the metallization of the 
Al(l) and Al(4) layers in actual samples is probably only 
partial, with a resulting decrease in tunneling width of a 
total of only one atomic layer. 



The LDOS for Co(4), which is the layer at the interface 
for the clean interface, shows the effect of hybridization 
with the s-p electrons in the Al(l) layer next to it. The 
d-electrons are more delocalized and this is reflected in 
that the d-electron peaks are not as well-defined as in 
the Co(l) layer. The peak in the spin- up d band near 3 
eV is reduced, and the peak near 5.6 eV has been shifted 
to a much smaller peak near 5 eV, and the rest of the 
spectral weight spectral weight has been moved to a more 
delocalized LDOS around 5 eV. 

There is also a slight reduction in the magnetic moment 
in this Co(4) layer as compared to the bulk-like Co(l). 



IV. INTERFACIAL SCATTERING AND 
FRIEDEL SCREENING EFFECTS 

From TEM studies of MTJ structures, it can be seen 
that there is atomic-level disorder at the interface ^ 6 ! 17 
This atomic level disorder at the interface will necessarily 
result in scattering of the conduction electrons. Assum- 
ing Al-termination at the interface, the Al atoms will act 
as positively-charged impurities, because of their charge 
in the alumina, and this charge should get screened by 
the s-p and d electrons in the Co. From Friedel's sum 
rule, we know that the conduction electrons will screen 
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Al(l)LDOS Plot 
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FIG. 4: LDOS for the first Al layer in the insulating barrier, 
Al(l), for the clean Al-terminated interface. 



the impurity charge, resulting in a change in the LDOS 
near the interface. Almost all of the screening will be 
due to the d-electrons as the LDOS of the d-electrons is 
much larger than the s- and p-electrons for Co at Ep. 
Also, since the spin-down DOS in cobalt is larger than 
the spin-up DOS at Ep , more spin-down electrons will 
be involved in the screening. This will result in a spin 
dependent change in the LDOS at the interface. Simi- 
lar spin dependent charge transfer effects have also been 
seen in other DFT works on clean Al-terminated inter- 
faces, where there is a larger transfer in the spin-down 
band than the spin-up band£ 

We can calculate the approximate shift in Ep and 
LDOS at Ep, Sp°(Ep), for the spin polarized d-electron 
bands, and give an estimate of the spin dependent 
Friedel-screening effect^ We assume one Al impurity 
per unit cell and use the LDOS of clean Co as shown 
in Fig. [2l The total phase shift for both spin-up and 
down bands must satisfy 9^(Ep) + 6i(Ep) = irZ, where 
Z is the charge of the impurity, so as to completely screen 
the impurity. The phase shift of each band, 9 s {e), is given 
by, 



tan((9 s (e)) 
F 8 {e) 



= P 



F s (e) - i 
P° s (e') 



de'. 



(2) 



Here, U, the impurity potential, as seen by both bands 
is clearly the same. Since p°(Ep) > p°(Ep), the spin- 
down band is expected to screen more of the impurity 
charge than the spin-up band. In the "rigid-band" ap- 
proximation, the shift in Ep can be approximated by 



A, = 



np°{E F ) ' 



and the change in the LDOS at Ep is 



approximated by 5p°(Ep) = A, 
ward calculation gives 9-[(Ep) 



^f(Ep). A straightfor- 
= £ 0.42tt and O^Ep) = 



0.58tt, which gives A T = 0.30 eV and A L = 0.56 eV, 
showing that the spin-down d-electrons experience a Frei- 
del screening effect that is about twice as strong. Since 
the tunneling current depends strongly on the DOS at 
Ep, this change in the LDOS will have a significant ef- 
fect on the spin polarized tunneling current. 

As will be shown below, DFT results for the Al-rich 
interfaces clearly demonstrates the significance of hy- 
bridization and the spin dependent Friedel screening ef- 
fect. In order to obtain an accurate and reliable band 
structures arising from this Al-induced disorder effect, we 
carried out DFT calculations on interfaces with an excess 
amount of Al incorporated at the interface. We believe 
the most obvious place for these excess Al atoms are to 
replace Co in the Co-lattice and we therefore created in- 
terfaces with Al-Co layers. As described in Section ITT1 we 
studied both junctions with one and two Co-Al interfacial 
layers. The case with only one interfacial layer can qual- 
itatively be interpreted as an interpolation between the 
clean and the two interfacial Co-Al layer interface so we 
will for clarity and compactness only report the results 
on the case with two interfacial Al-Co layers. The layer 
spin polarization in Fig. I1Q[ which includes the results 
for all three structures, clearly shows the intermediate 
features for the one interfacial layer structure. 




FIG. 5: Atomic structure of a disordered Al-rich interface 
junction. The unit vectors for the supercell as well as the 
layer identification tags are shown. The Int25 layers contain 
a 3:1 Co to Al ratio whereas the Int50 layers have a 1:1 ratio. 
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V. AL-RICH DISORDERED INTERFACE 

The self-consistent atomic structure of the supercell 
for modeling the Al-rich disordered interface, with inter- 
face layers Int25 and Int50, is shown in Fig. [5] The first 
interface layer Int25(l), i.e. the layer closest to the Co 
substrate, has a 3:1 ratio of Co to Al, and the second 
mixed layer Int50(l), which is next to the Al(l) layer, 
has a 1:1 ratio of Co to Al. In Figs. EE] we plot the 
spin dependent LDOS for the first Co layer, Co(l), the 
first interface layer, Int25(l), the second interface layer, 
Int(50), and the first Al layer, Al(l), respectively, to show 
the change in spin dependent LDOS as we approach and 
cross the interface. Comparing these results with those 
from Sec. IIII1 we are able to see the effects of Friedel 
screening of Al atoms at the disordered interface. 



Interface Layer 1 (Ca^AL.J LDOS Plot 
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FIG. 6: Spin-dependent LDOS for Co(l) layer. 

As expected, the LDOS for Co(l) is almost the same 
as bulk Co and the Co(l) layer of the clean interface case, 
giving an exchange splitting of w 1.5 eV. The LDOS of 
Co(l) shows a large spin-down polarization at Ep 7 with 
Pl(Ep) ~ 1 eV -1 and pi(Ef) ~ 4 eV -1 , similar to that 
of the clean interface case. This means the screening 
length for the intcrfacial disorder effects is very short, of 
about one atomic layer, as expected for a metal. 

The mixing of Al and Co in the disordered layers near 
the interface results in spin dependent Friedel screening, 
which changes the LDOS. Since pi(Ep) > p^(Ep) for 
bulk Co and also in the Co(l) layer, screening of the 
Al atoms in the Int25 and Int50 layers would be mostly 
from the spin-down d-electrons. Therefore pi(Ep) should 
change proportionally more than p^(Ep), as shown in the 
calculation for the Friedel screening. 

The effects of the spin dependent screening is seen in 
both the LDOS near Ep and the total number of d elec- 
trons per layer. There is a large change in the LDOS at 
Ep for the Int50 layer compared to the Co(l) layer, which 
has significant importance for tunneling transport. For 



FIG. 7: Spin-dependent LDOS for first interface layer, 
Int25(l), with the ratio Co:Al = 3:1. 

Interface Layer 2 (Co,„„AL.„) LDOS Plot 
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FIG. 8: Spin-dependent LDOS for second interface layer, 
Int50(l), with the ratio Co:Al = 1:1. 



the Int50 layer, Pd\ — 1.0 eV x , and 



1.9 eV" 



Comparing this to the bulk- like layer of Co(l) where 
Pd\ = 1.0 eV -1 and p d ± — 3.4 eV -1 , we see that there 
is clearly a much larger change in the spin-down LDOS. 
The spin-up LDOS is unchanged, while the spin-down 
LDOS is almost reduced by half. Similarly, there is a 
much larger change for the total spin-down d electrons 
in the Int50 layer. The total d electron charge in the 
Int50 layer is 8.4 eV _1 /layer and 7.2 eV _1 /layer for 
spin-up and down respectively. In comparison, the av- 
erage of the Al(l) and Co(l) layers is 9.0 eV _1 /layer and 
6.0 eV _1 /layer respectively for spin-up and down. This 
shows the spin dependent change in the LDOS of the 
interfacial layers is due to screening effects, whereas hy- 
bridization should lead to a similar change in the LDOS 
for both spins. 
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In addition, we have calculated the average of the 
LDOS of the Co(l) and Al(l) layers in the clean case 
to obtain an estimate for the effects of hybridization, 
and compared that to the Int50 LDOS for a clearer pic- 
ture of the effects of the spin-dependent screening. The 
LDOS near Ef for the two are very similar for spin-up, 
whereas spin-down shows a large shift of spectral weight 
from around Ep to the peak seen in Int50 at about 4.6 
eV. This leads to the change in the spin polarization at 
Ep and also affects the spin polarized tunneling current, 
both of which are discussed further in the following sec- 
tions. 

The sharp peaks in the d-electron LDOS that were 
seen in the Co(l) layer are also significantly reduced in 
the Int25 and Int50 layers, indicating that the d-electron 
orbitals are more delocalized. This is due to hybridiza- 
tion with the Al s-p bands from the neighbouring Al(l) 
layer, and with the Al s-p electrons in the layers itself. 
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FIG. 9: Spin-dependent LDOS for first Al layer, Al(l). 

We notice that the first Al layer, Al(l), has a small 
but finite LDOS at Ef, since it has become metallized 
by proximity to the Co layers, an effect also seen for 
the clean interface. The effective insulating width of the 
AI2O3 layer is therefore again reduced by approximately 
one to two atomic layers. 

We plot in Fig.rTU]the spin polarization, S s , of the three 
different Al-terminated MTJ structure studied, starting 
in the bulk and moving across the interface into the 
AI2O3 insulating barrier. As we approach the interface 
we can see an increase in S s of each layer. The bulk Co 
layers have a, 6 S w —(0.6 — 0.8) in all structures. For 
both the clean and disordered system, 5 S increases as we 
approach the interface. For the clean Al-terminated in- 
terface the spin polarization remains slightly negative at 
the interfacial Al(l) layer, and becomes slightly positive 
in the 0(1) layer. For the two interfacial layer struc- 
ture (Fig. [5]) the spin polarization similarly increases to- 
wards the interface, 6 S = —0.32 in the Int25 layer, and 
5 S = —0.27 in the Int50 layer. Furthermore, the spin po- 
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FIG. 10: Spin polarization at Ep, 8 a (Ep), of the three dif- 
ferent MTJ structures investigated assigning the Al(l) layer 
to be layer 0. Clean interface (black), one Co/Al interfacial 
layer (green), and two Co/Al interfacial layers (red). 



larization becomes positive, S s = 0.2 , at the metallized 
interfacial Al(l) layer. 

Looking at the change in spin polarization from the 
bulk- like Co(l) layer to the layer at the interface, the in- 
crease in 6 S from the Co(l) to Co (4) layer for the clean 
interface system can be attributed solely to hybridiza- 
tion with the Al s-p electron in the Al(l) layer. For this 
system, the spin polarization changes from -0.64 to -0.43, 
which is a change of 33%. On the other hand, for the dis- 
ordered system with two mixed layers, S s increases from 
-0.58 to -0.26 between the Co(l) and the Int50 layer, 
which is a change of 55%. This clearly shows the addi- 
tional effects of the spin dependent Friedel screening and 
interfacial disorder on the spin polarization. 

To study the source of the positive spin polarization 
in the interface region shown in Fig. I10[ we have per- 
formed a calculation on a model system, an 'essential 
unit cell' of the interface region, shown in Fig. [TTJ This 
is composed of an Al atom bonded to two Co atoms in 
one direction, and to three oxygen in the other. These 
oxygen are bonded to a second Al, and the structure 
simplified by termination at this point with two final Co. 
The bond lengths and angles were chosen to match those 
of the clean structure, Fig. LTJ For reasons shown in the 
following discussion, we believe the model structure may 
capture much of the essence of the spin interactions and 
hybridization at the interface. 

Density functional calculations were performed on this 
system in a comparable manner to those in the rest of 
this paper, including full relaxation of the structure. The 
program Quantum Espresso^ was used, with conver- 
gence and other properties similar to those described in 
Sec. HIl Relaxation of the structure did not change the 
bond lengths or angles from those of the full calculation 

of Fig. ru 
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terminated systems shown in Fig.fTOlshows explicitly that 
not only a O-rich but also a Al-rich disordered interface 
can give rise to a positive spin polarization in the insu- 
lating alumina layer. Moreover, this positive spin polar- 
ization can be reached in a region with a finite LDOS 
at the Fermi level for an Al-rich interface, namely the 
Al(l) layer, which is of crucial importance when study- 
ing tunneling processes. The effects of the change in spin 
polarization on the tunneling current is modeled using 
spin dependent effective masses for the different bands in 
each layer as discussed the section below. This is shown 
in Table. HI especially for the d2 bands, which contribute 
the most to the LDOS at Ep. The details of the effective 
mass derivations, and effects on the tunneling current are 
discussed in the following Sec. IVfl 




FIG. 11: Atomic structure of model system of 9 atoms repre- 
senting the Al-terminated C0/AI2O3 MTJ structure. The fig- 
ure shows the total spin polarization, and the color indicates 
the spin polarization for each atom. This result, matching 
so closely the larger calculation, is suggestive that superex- 
change is the mechanism for the positive spin polarization in 
the larger structure, and is definitely the mechanism for the 
observed spin polarization in this smaller structure. 



VI. SPIN POLARIZED TUNNELING CURRENT 

In the section above we showed that Al-terminated and 
Al-rich interfaces can give rise to a net positive spin po- 
larization in the insulating barrier. We now explicitly 
show that the tunneling current is also positively spin 
polarized. In order to calculate the tunneling current, we 
model the MTJ system by a position-dependent effective 
mass model, which we then solve exactly to obtain the 
tunneling current. The tunneling layer is modeled in the 
standard manner as a potential barrier, V , with thick- 
ness a, and applied voltage V\ . The Hamiltonian for the 
system is written simply as an effective kinetic mass term 
and a potential barrier: 



Fig. ITTI shows the spatial distribution of the spin polar- 
ization, with contrasting color indicating opposite direc- 
tion of spins. We note first of all the marked similarities 
between these results and those of Fig. [TUJ The total 
spin polarization on both the Al and O layers is positive, 
opposite to the negative of the Co. In addition, the spin 
polarization of the Al is not centered on the Al itself, 
but rather is midway to the Co, indicating the extent 
of the hybridization. Both the Al and Co polarizations 
have a marked anisotropy, indicating the preferential in- 
volvement of certain orbitals in the development of the 
spin polarization. Al draws electrons from the Co, a pro- 
cess akin to Friedel screening in a larger system. The 
electrons near the Al are spin polarized in the opposite 
direction to Co due to superexchange in this small sys- 
tem. Similarly we find in the structure in Fig. [TJ that 
the Al has a total spin polarization opposite to the Co 
atom. We suggest this is due to similar Pauli exclusion 
effects in the screening process. If the MTJ structure is 
O-terminated instead, the interfacial O atoms would be- 
have similarly to the interfacial Al-atoms that we have 
described above; i.e. the O atom would be screened by 
the Co electrons and would have a net positive total spin 
polarization. 

The spin polarization results for the three different Al- 



H = + 6(x)6(a - x)(V - —x). 

Zm(x) a 



(3) 



The effective mass, m(x), changes discontinuously 
from layer to layer, and the values for the different s-p 
and d bands are obtained from the DFT calculations as 
described below. The Fermi energy is 6.78 eV, as deter- 
mined from the LDA results, the effective barrier height, 
Vq = 9.4 eV, is obtained from electron spectroscopy of 
thin film A1 2 3 layers^ The effective width of the A1 2 3 
layer, a, is taken to be one atomic layer thinner than the 
actual thickness shown in Fig. due to the partial metal- 
lization of the Al layers at the interface. The appropriate 
boundary conditions to be imposed between every layer 



is, 



V(o-) = v(o+) 



1 d ^ (0 
mi dx 



1 

777,2 dx ' 



(4) 



where 771 1 is the effective mass for the layer on the left, 
x < 0, and 7772 is the effective mass for the layer on the 
right, x > 0, where x = is where the layers meet. 

We extract effective masses for our model from the 
DFT results as follows. The DFT LDOS data is pro- 
jected into s-p and d bands for each of the layers, which 
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we model separately with different effective mass values. 
The d-electron band is modeled using three bands, dl, d2, 
and d3, so as to more realistically capture the LDOS near 
Ep. First, we project the DFT LDOS data into s-p and 
d bands for each of the layers and in particular look at 
the data near Ep. Next, we model the d-electron band as 
three separate bands in order to more accurately capture 
the LDOS near Ep. The large peak near Ep is modeled 
using one band, labelled as d3, the overall large-effective 
mass d band is modeled using another band, labelled as 
d2, and a more flat sp-like band is labelled as dl. We 
model the s-p, dl, d2, and d3 bands with a free-electron 
density-of-states expression, with effective mass, m*, and 
the band edge as fitting parameters. 

To check our results, we also use the full LDA band 
structure to extract the Fermi velocity, vf, for each band 
and then use Vf = to extract an effective mass, m* for 
each of the bands. We find that the effective masses, m, 
extracted from the free electron picture above are in good 
agreement with the band structure masses, m* , and thus 
we have a high confidence in our approach. The effective 
masses, with respect to the bare electron mass mo, are 
shown in Table |TJ Here the two disordered layers at the 
interface, Int25 and Int50, are treated as one single layer 
for ease of calculation. The LDOS of the two layers are 
similar, and the effective mass value used is the average 
of the effective masses obtained from the two layers. The 
d3 bands in the interface layers do not cross Ep and 
are therefore assumed to be irrelevant for the transport 
calculations. 





Co(l) Layer 


Int25/Int50 Layer 


AI2O3 Layer 


Bands 


s-p 


dl 


d2 


d3 


s-p 


dl 


d2 


s-p 


Spin-t 


0.45 


0.81 


2.05 


6.19 


0.51 


0.81 


1.70 


1.0 


Spin-| 


0.37 


1.02 


4.99 


7.89 


0.51 


0.58 


1.70 


1.0 



TABLE I: Effective mass values (m/mo) for the different elec- 
tron bands in the cobalt layer, the interface layer, and the in- 
sulating AI2O3 layer. Notice the large change in the effective 
mass for the spin-down d2 bands between the Co(l) and the 
Int25/Int50 layers. 



The increase in the spin polarization due to the Al- 
rich interface is reflected in the spin dependent effective 
masses of the d-electron bands, especially in the d2 bands 
which has the largest LDOS at Ep. It is shown in Ta- 
ble U that the d2 bands in the Co(l) layer have effective 
masses 2.05 mo and 4.99 mo for spin-up and down respec- 
tively, which both reduce to 1.70 mo in the interface layer. 
Clearly the larger decrease in spin-down LDOS near Ep, 
due to the Friedel screening effect and LDA results shown 
in Sees. [TV] and El leads to a corresponding increase in 
spin polarization 5 s (Ep). This decrease in LDOS is mod- 
eled by the change of the effective mass values from the 
Co(l) layer to the interface layers. A smaller tunneling 
probability for the spin-down d bands hence results, due 
to the larger spin-down effective masses, and the greater 
mis-match in effective masses between the Co(l) and the 
Int25/Int50 layer. This is clearly seen for all three dl, 



d2 and d3 bands, and hence there is a smaller spin-down 
polarized tunneling current. 

The exact wavefunctions for the Hamiltonian, Eq. ([3]), 
are described by Airy functions whose coefficients are 
determined numerically using the boundary conditions, 
Eqs. (01), and a normalization condition for the wave- 
function. The numerically determined wavefunction is 
then used to calculate the tunneling probabilities, and 
tunneling current which is given by the standard current 
expression, 

j=in (?m r{SB) -^^ x) y ( 5 ) 

Figure [T2l shows the tunneling current for the parallel 
configuration, where the magnetization of the two cobalt 
layers are parallel. Therefore, the spin-up electrons tun- 
nel to the spin-up bands, and similarly for the spin-down 
electrons. The s-p band electrons naturally tunnel to 
the s-p band by symmetry, and since the dl, d2 and d3 
band electrons were originally modeled from the same d- 
electron band, they tunnel into the other bands as well; 
i.e. the dl band electrons tunnel into the d2 and d3 bands 
as well as the dl band. 

Cobalt Spin Polarized Current 



s & d Bands Contribution to Spin Polarized Current 




Voltage (V) 



FIG. 12: Spin polarized tunneling current of MTJ structure. 
The lower transmission coefficient for the spin-down bands, 
due to a greater mis-match in effective masses, result in a 
lower tunneling current compared to the spin-up bands. 

The spin polarized tunneling for both the s-p electrons 
and the d-electrons are positive, up to s» 0.9 eV. The 
positive spin polarization tunneling current arises from 
a greater mis-match in the effective masses of the spin- 
down bands between the bulk cobalt layer, Co(l), and 
the interface layers as compared to a smaller mis-match 
in the spin-up bands. This is also seen in the s-p bands, 
where there is a slightly better match between the spin- 
up effective masses. This results in a slightly larger spin- 
up tunneling current for the spin-up s-p band, as shown 
in Fig. [TJ] There is also a large peak, modeled as the d3 
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band, in the LDOS of spin-up d-electrons, approximately 
100 meV below Ep, that affects the spin polarized tunnel- 
ing; whereas a similar peak for the spin-down d-electrons 
is 1.26 eV above Ep due to the exchange splitting. This 
is seen clearly in Fig. [6] Therefore, for an applied volt- 
age larger than 100 meV, this large spin-up d3 band will 
contribute, and results in a larger spin-up polarized tun- 
neling current. This continues until 0.91 eV, when the 
bottom of the spin-down d3 band crosses the quasi- Fermi 
level and starts to contribute to tunneling. The positive 
spin polarized tunneling current reflects the positive spin 
polarization shown in Fig. [TO] as we cross from the cobalt 
layers into the AI2O3 layer, and explains the positive spin 
polarization found experimentally. 

Finally, we can numerically calculate the tunneling 
spin polarization, P = j—^, where J-j-j is the spin polar- 
ized tunneling current given in Fig. [TOJ Using Julliere's 
formula, Eq. (Q]), we then calculated the TMR from the 
numerically derived tunneling spin polarization P, and 
the results are shown in Fig. [13] 

Tunneling Magnetoresistance 

TMR - (R Ap - RpVfRy.) 
0.5 1 , 1 1 1 , 1 , 1 1 1 , 1 , 1 1 1 , 1 , I 




Voltage (V) 

FIG. 13: Tunneling magnetoresistance of MTJ structure cal- 
culated from Julliere's formula, with spin polarization ob- 
tained from spin polarized tunneling current shown in Fig. 1121 

The calculated spin polarization is P — 53%, which 
compares very well with the experimental values of 42% 
for Co and 55% for Co 50 Fe 5 ci£2i The calculated TMR is 
43% at zero applied voltage, which compares favourably 
with the experimental value of 30%. Next, we see that 
the TMR decreases monotonically with applied voltage 
and reaches zero at 0.91 V. The spin polarization and 
TMR both decrease to zero at this voltage, which is the 
spin-down d3 band edge and is of the order of the ex- 
change splitting. Finally, we also see two well-defined 
"roll-overs" in the TMR curve labeled "d2| band edge" 
and u d2± band edge". This is due to band structure 
effects; the first "roll-over" at 0.23 V is due to the spin- 
up d2 band edge and the second "roll-over" at 0.375 V 
is due to the spin-down d2 band edge, where the two 
bands stop contributing to the tunneling transport. The 



monotonic decrease of TMR with increase in voltage and 
the "roll-overs" have been seen experimentally, see for 
e.g. Refs.iiHi These two "roll-overs" are much more 
prominent in the numerical results because the band 
edges for the c?2| and d2[ bands in the effective-mass 
model are modeled using step-functions, which is more 
abrupt than the real band structure. 



VII. CONCLUSIONS 

In this paper, we have studied the effect of interface 
mixing and atomic-level disorder on the spin polarized 
tunneling current in C0/AI2O3 MTJs. The spin polar- 
ization of tunneling currents in MTJs have clearly been 
seen in experiments to be positive, at odds with simplistic 
tunneling models which depends on the spin polarization 
of the DOS at the Fermi energy of the leads. There is 
now a growing body of evidence that interface effects at 
the boundary between the Co layers and the aluminum 
oxide insulating layer play an important role in determin- 
ing the spin polarized transport in these systems. A large 
number of studies have assumed an O-rich interface, and 
found positive spin polarization for such systems. How- 
ever, we believe that an Al- rich/terminated interface is 
more relevant experimentally when the junction is not 
heavily oxidized. 

The MTJ was modeled realistically using DFT in both 
the LDA and GGA approximations, of which the LDA 
scheme was found to be more consistently reliable. We 
have studied both clean Al-terminated and Al-rich dis- 
ordered interface systems, and found a consistent trend 
of increasing positive spin polarization with Al content 
of the interface. The main effects of the Al-rich disor- 
dered interface are spin dependent Freidel screening of 
the Al-impurities in the interface lead layers, and hy- 
bridization with the s-p electrons of the Al-impuritics 
and the neighbouring first Al layer of the insulating bar- 
rier. This causes a much larger change in the spin-down 
d-electron LDOS near Ep of the interface layers, leading 
to an increase in the spin polarization near Ep. An- 
other important result of our LDA calculations is that 
the proximity effect leads to metallization of the first Al- 
layer by the neighbouring Co layer. This reduces the 
effective width of the insulating layer, and can be seen in 
the slope of the TV curve. 

Studying the interface with a model system which fo- 
cuses on the Co-Al-0 interactions, we find good agree- 
ment with the larger calculation on the incidence and sign 
of spin polarization in the interface layers. This model 
calculation suggests that the source of the opposite po- 
larization of the Al and O is a superexchange mechanism 
similar to that seen with magnetic atoms on an insulat- 
ing surfaced Just as the RKKY interaction was origi- 
nally studied with isolated magnetic impurities in bulk 
metals, and then widely seen in metallic magnetic multi- 
layers, we postulate that the superexchange interactions 
previously seen for isolated magnetic atoms separated by 
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a nonmagnetic insulating 'island' are also applicable to 
many magnetic tunnel junctions. Because the hopping 
interaction dies down with distance into the insulator, 
the net effect will be an opposite spin polarization of the 
interface layers. 

The effects of increasing spin polarization in the LDOS 
near Ep is captured using a layer and spin dependent 
effective mass model, which is used to calculate the tun- 
neling current. The greater mis-match of the effective 
masses leads to a smaller transmission coefficient for the 
spin-down d bands, and the tunneling current clearly 
shows a positive spin polarization that is in agreement 
with experimental results. This defines an effective tun- 
neling spin polarization, which allowed us to calculate 
the TMR using the Julliere model, giving results that 
are also in good agreement with experiments. 

In summary, using DFT calculations, we have found 
that Al-terminated interfaces in MTJs show positive spin 
polarization of the LDOS at Ep, showing that positive 
spin polarization is possible not only in O-terminated 



MTJ systems. Furthermore, we formulated a spin depen- 
dent effective mass model, using values extracted from 
the DFT results, that was used to calculate the tunneling 
current and the TMR. The results clearly show a positive 
spin polarization and TMR values that are in agreement 
with experiments. Thus, our model demonstrates how 
positive spin polarization and realistic TMR values can 
be obtained in an Al-terminated MTJ. 
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